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Abstract. The Bose- Hubbard Hamiltonian capturing the essential physics of the 
arrays of interacting Bose-Einstein condensates is addressed, focusing on arrays 
consisting of two (dimer) and three (trimer) sites. In the former case, some results 
concerning the persistence of mean-field features in the energy spectrum of the 
symmetric dimer are extended to the asymmetric version of the system, where the two 
sites are characterized by different on-site energies. Based on a previous systematic 
study of the mean-field limit of the trimer, where the dynamics is exhaustively described 
in terms of its fixed points for every choice of the significant parameters, an interesting 
mapping between the dimer and the trimer is emphasized and used as a guide in 
investigating the persistence of mean-field features in the rather complex energy 
spectrum of the trimer. These results form the basis for the systematic investigation of 
the purely quantum trimer extending and completing the existing mean-field analysis. 
In this respect we recall that, similar to larger arrays, the trimer is characterized 
by a non-integrable mean-field dynamics featuring chaotic trajectories. Hence, the 
correspondence between mean-field fixed points and quantum energy levels emphasized 
in the present work may provide a key to investigate the quantum counterpart of 
classical instability. 



PACS numbers: 03.75.Kk, 03.65. Sq, 05.45.Mt 
1. Introduction 

In a recent publication pQ , the present authors reported an exhaustive study of the mean- 
field (MF) dynamics of a system of three interacting Bose-Einstein condensates. Such 
a system, henceforth referred to as trimer, deserves attention since it can be considered 
a paradigm for larger arrays. Indeed, similar to longer chains and more in general to 
larger arrays, its MF dynamics is non-integrable and its phase space is characterized by 
extended regions of instability On the other hand, the system is still sufficiently small 
to allow a thorough analytical study, which provides a valuable operational tool in both 
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interpreting experimental observations and identifying significant experimental setups 
triggering macroscopic effects PQI21- Furthermore, such exhaustive investigation makes 
the trimer a candidate benchmark for verifying the effectiveness of the MF description, 
which, being essentially equivalent to the Bogolubov approximation, is expected to 
provide a reliable picture of the dynamics of an array of interacting BECs for sufficiently 
large boson populations. 

As to the experimental realization of interacting BEC arrays, optical fragmentation 
allowed to attain ID 01 E], 2D © and even 3D [2j lattices consisting of a rather 
large number of sites. Though Bose-Einstein condensation has been already achieved 
both in an optomagnetic jH] and, more recently, in an entirely magnetic [3 QDj double- 
well trap, microtrap technology ^IJ ^] is probably the most promising approach for 
realizing the small systems addressed by the present article, as well as any complex array 
characterized by a trapping potential lacking the typical periodicity of optical lattices. 

We will focus on two small systems, the asymmetric dimer, i.e. a condensate 
trapped into an asymmetric double well trap, and the open trimer of reference £Q, where 
an energy offset between the central and the lateral wells of the trapping potential is 
present (a multiwell trapping potential characterized by uneven well depths may possibly 
be achieved within a supercell of a suitable optical superlattice [EH El)- Our aim is to 
compare the spectrum of the quantum Hamiltonians describing these systems with the 
characteristic energies of the relevant mean- field picture, in the spirit of reference [TKj . 
There Aubry and co-workers illustrated a remarkable relation between the energy of the 
mean-field saddle point (separatrix) of the symmetric dimer and some features of the 
spectrum of the quantum counterpart of the same system. 

In the following we will extend these results to the more general case of the 
asymmetric dimer. Subsequently, exploiting the exact mapping between the mean-field 
picture of the dimer and that of the integrable sub regime of the trimer, which we show 
to apply at a qualitative level in the purely quantum treatment of the systems, we will 
point out some more interesting manifestations of MF features in the spectrum of the 
trimer. 

Before proceeding with our analysis, we briefly recall some generalities. The 
essential physics of an array of coupled BECs is captured by the Bose-Hubbard model 
|16j . describing a gas of interacting bosons hopping across the sites of an ambient lattice. 
The relevant Hamiltonian reads H = J2j[Un* — VjUj — T/2 J2k~j( a t ak ~^ a fc a j h )]> where 
off (a,j) is the creation (annihilation) operator on site j, nj = a^aj counts the number 
of bosons at the same on site, U is the on-site interaction parameter, Vj is a the energy 
offset of site j, T is the hopping parameter and the subscript k ~ j restricts the sum to 
the nearest neighbours of site j. 

A mean-field treatment of the Bose-Hubbard model can be achieved by means of a 
time dependent variational principle based on Glauber's coherent states \Tj\. According 
to such picture each site in the array is described by a macroscopic complex variable 
Zj = ^fp~j^ i whose square modulus pj and argument <pj represent the population 
and phase associated to the condensate sitting at that site, respectively. Barring 
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constant terms, the effective MF Hamiltonian has the form 7i = YljW\ z j\ ~~ v j\ z j\ ~ 
T/2 J2k~j( z j z k + z j z k)) anc ^' equipped with the Poisson brackets among the (canonical) 
macroscopic variables, {zj,zl} = 5jk determines the dynamics of the system. ||. The 
essential features of the dynamics are determined by the phase space structure, which 
can be analyzed working out the location and the stability character of the corresponding 
fixed points. Notice that since H features only two-boson operators, it commutes with 
the the total number operator n = £\ nj. Likewise, the total population p = £\ pj is a 
constant of the motion mirroring the global phase invariance of TC. Hence the fixed point 
equations have the form {zj, 7i — fi p} = 0, where fi is a Lagrange multiplier constraining 
the total number of bosons. Though the solutions of the fixed point equations are in 
principle complex vectors, they can be safely restricted to real values, with the only 
exception of the topological solutions [HUE] which are to be found in cyclically invariant 
structures (i.e. homogeneous closed chains and their higher dimensional generalizations). 
Since the structures we will consider are neither cyclic nor homogeneous, in solving the 
fixed point equations we will set Zj = y/pXj, where Xj £ R and J2j x | = 1- 

2. The asymmetric dimer 

The Bose-Hubbard Hamiltonian for an asymmetric dimer has the form 

TV T 

H d = U d (n\ + n 2 r ) - — (m - n r ) - — {at a+ + aj a r ) (1) 

where the subscripts £ and r denote the left and the right sites of the structure, 
respectively, and w d is a local energy offset deriving from the asymmetry of the double- 
well trap confining the BEC. As we briefly recalled above, the MF counterpart of 
Hamiltonian © is H d = U d (\z e \ 4 + \z r \ 4 ) - ±w d (\z e \ 2 - \z r \ 2 ) - \T d (z lZ * r + z* e z r ), 
and determines the evolution of the macroscopic semiclassical variables through the 
(integrable) equations of motion 

itet = (2U d \z e \ 2 - z e - y z t, iter = (2U d \z r \ 2 + z r - y z £ (2) 

According to the discussion in the previous section, and after introducing the 
conveniently rescaled parameters 

w d T d Xd ( o\ 

= 77 — , r d = - — , /x d = - — (3) 
U d p d U d p d U d p d 

where p d = \zi\ 2 + \z T \ 2 = p d (xj + xf), the consequent fixed point equations read 

2x 1 ~ y ~ A*d) x t - yx r = 0, (2x 2 r + y - /i d ) X r - yX^ = (4) 

Notice that r d and u d are the only significant parameters, whereas /u, d has an auxiliary 
function and its value is different for different fixed points. In the following we will 

|| The set of equations governing the dynamics of the system are basically the same discrete nonlinear 
Schrodinger equations (DNLEs) deriving from the application of the tight-binding approximation to the 
Gross-Pitaevski equation describing the evolution of a fragmented BEC ^21 • Reference reports on 
some analytical solutions of similar DNLSEs on small arrays such as the ones presently under concern. 
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consider |z/ d |'s up to the order of 2, corresponding to relative energy offsets of the order 
of 2 U pa- If large total populations are involved, the limit situation may correspond to 
an overly asymmetric double well potential. Notice however that comparable relative 
energies are easily attained in the MF integrable subregime of the trimer, which can be 
described in terms of an asymmetric dimer, as we will discuss in the next section. 

Depending on whether \ud\ is smaller or larger than (2 2 / 3 — t^ 3 ) 3 / 2 , Q have either 
four or two solutions. In the former case there are two maxima, a saddle point and a 
minimum. As |fd| grows, the local maximum and the saddle get closer and closer, until 
they collapse and disappear when |z/ d | crosses the (r d -dependent) critical value ^f. We 
mention that the presence of two maxima is an aspect of the mirror-symmetry breaking 
induced by nonlinearity, leading to interesting phenomena such as self-trapping and 
7r-phase oscillations, that are investigated in reference j2l] for the symmetric dimer. 

Aubry and co-workers [THj focused on the symmetric dimer (u>d = =^> v& = 0), 
and found that when the MF saddle point (separatrix) is present, i.e. for r d < 2, its 
energy is strictly related to some interesting features of the spectrum of Hd- They noticed 
that the density of states of H d features a sharp peak centered around the saddle point 
energy. Furthermore this MF energy marks the boundary between two regions of the 
spectrum characterized by significantly different magnitudes of the splitting energy of 
the doublets arising from the mirror symmetry of the symmetric dimer. More precisely 
the energy splitting experiences a drop of several orders of magnitude as the average 
energy of the doublet becomes larger than the threshold value marked by the saddle point 
energy. Further analytical considerations about the splitting pattern of the quantum 
dimer doublets and its relation with the saddle point energy of the corresponding mean- 
field model can be found in reference [22]. Notice that the magnitude of the splitting 
between even-odd doublets of the symmetric dimer is also related to the characteristic 
time-scale for the symmetry-restoring quantum tunneling between semiclassically self- 
trapped states p5 | l26 t l22]. 

Here we show that the results by Aubry and co-workers can be extended and 
generalized to the asymmetric dimer, where Wd (and of course z/d) has a finite value. 
To this end it proves useful to plot the energies of the MF fixed points of H d versus 
the energy offset among the two sites of the dimer, Wd, which gives a measure of its 
asymmetry. This can be accomplished recasting the dynamics of the trimer in terms 
of the alternative set of canonical variables (See e.g. reference [27] for further details 
about this alternative canonical representation), pd = \zi\ 2 + \z2\ 2 , D = \zi\ 2 — \z2\ 2 , 
9 = (0i — </>2)/2, ip = (</>i + 02)/2, whose significant Poisson brackets are {D,6} = 
{pd,4>} = ft -1 - This way, the relevant dynamic equations can be used to give a 
parametric representation allowing to plot the fixed point energies versus Wd for a 
given value of the hopping amplitude. Since it is interesting to compare situations 

% Notice that the coalescence of the (possibly local) maximum and the saddle can be likewise achieved 
also through a variation of rd, and occurs at a | ^-dependent critical value ^J|52]. A similar coupling- 
driven change in the semiclassical phase space of a 2D BEC array has been recently proposed as a 
mechanism to trigger an abrupt derealization in the system |23| . 
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Figure 1. Comparison between the mean- field and quantum characteristic energies 
for the asymmetric dimer for = 0.1. Left panel. Thin solid lines: energy levels 
for a dimer containing pd = 16 bosons. Alternating colours were used to emphasize 
the avoided level crossing. Thick black line: MF minimum; Thick black dashed lines: 
MF maxima. Circles: MF saddle. Right panel. The MF energies are compared to 
the density of the states for a dimer containing p& = 200 bosons. The inset shows the 
spectrum and the energy of the MF saddle in a small region of the figure plane. 



characterized by different total populations, we provide this parametric representation 
in terms of the rescaled parameters Td, u d , 

. , / r d cos2#\ . (l-d 2 )l -r d cos2# 

where E d is the energy in units of Up d , d = D/p d e [—1, 1] is the rescaled population 
imbalance [23|, and 9 is either zero or tt/2. The former choice pertains to the minimum 
fixed point no matter for the value of the hopping amplitude. As we discussed above, 
when Td > 2 there is only one more point, a maximum, independent of the value of 
the energy offset \vd\- m this case © the remaining choice for the angular variable, 
9 = 7r/2, describes the energy of the maximum on the entire domain, d G [—1, 1]. If, 
conversely, Td < 2, the choice 9 = n/2 pertains to different fixed points depending on 
the domain of evaluation. More in detail, the energy of the saddle, local maximum and 
absolute maximum are obtained evaluating (JHJ) with \d\ belonging to the interval [0, d±], 
[di, d 2 ] and [d 2 , 1], respectively, where d x = y/l — (r d /2) 2 / 3 and d 2 = a/1 — (r d /2) 2 . 

The circles, the thick solid line and the thick dashed lines appearing in both panels 
of figure d were obtained as described above, and show the behaviour of the energy 
of the MF saddle, minimum and maxima, respectively, for Td = 0.1. In the left panel 
of the same figure we also plotted the energy levels of H d for a total population of 
Pd — 16 bosons (thin solid lines), using alternating colors for emphasizing the avoided 
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level crossing feature. Notice indeed that at a discrete set of values of the energy 
offset (wd = 2kUd, k — 1, . . . ,Pd), signalled by the vertical dotted lines, some pairs of 
adjacent energy levels become so close that, if plotted with the same color, they would 
give the illusion of crossing as they actually do when r d = 0. If conversely r d > 0, 
the adjacent energy levels never cross" 1- , but rather they "bounce" on each other, so that 
their behaviour with varying v& resembles the trajectory of particles repelling each other 
(see e.g. [2H])- Hence, whenever the energy offset between the wells is an even multiple 
of Ud (^d = 2/c/pa), the situation becomes similar to that of the symmetric dimer: the 
couples of degenerate levels possibly present for r d = split for finite values of the 
hopping parameter. Once again, the extent of the splitting of a given doublet can vary 
of order of magnitudes depending on whether its average energy is smaller or larger 
than the MF saddle point energy. Hence also in the asymmetric case such energy marks 
a clear boundary between two regions of the spectrum with quite different features. 
This is rather evident in the inset of the right panel of figure Q where we plotted the 
energy levels for a quantum dimer with a total population of p d = 200 bosons (black 
and gray solid lines) in a small region around the saddle point energy (white circles). 
The main plot of the right panel of figure Q shows the (Lorentzian broadened) density 
of states of Hamiltonian (JTJ for the same total population. It appears evident that also 
in the asymmetric case the sharp peak of the density of states (dark shading) follows 
very closely the saddle point energy (white filled circles). These features are observed 
whenever the hopping amplitude is sufficiently small (r d < 2, T d < 2[/ d p d ) to allow the 
appearance of a saddle point for some values of Vd- 

3. The trimer 

In reference pQ the present authors analyzed the mean-field dynamics of a trimer 
described by the following Bose-Hubbard Hamiltonian 

3 ^ rj-\ 

Ht = U t J2 n ) + y ( n i ~ n 2 + n 3 ) - — [(ax + a 3 ) + (af + a£) a 2 ] (6) 

where both the relative offset of the central well, w t , and the hopping amplitude T t 
were regarded as adjustable parameters. The results of this study were gathered into 
two stability diagrams summarizing, for any choice of the significant parameters, the 
character and linear stability of the entire set of MF fixed points, and hence describing 
the essential features of the phase space structure. This knowledge would allow to 
produce the trimeric counterparts of figure ^ with no further effort other than the 
diagonalization of Hamiltonian ©. Notice however that the size of this operator 
scales as (p t + l)(pt + 2)/2, which could be so large to make analytical diagonalization 
impracticable even for a relatively small total boson population p t . Thus the only readily 

+ We mention that the avoided level crossing is discussed to some extent in reference |28| . However 
there the varying parameter is the hopping amplitude between the sites of a dimer and a third site 
lacking of the repulsive term. 
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available analytical information to provide a systematic description of the quite intricate 
spectrum of the trimer (see e.g. figure EJ) is the Z 2 symmetry of (jUJ), allowing to classify 
its energy levels according to their parity. 

This is why here we prefer to proceed discussing a close relation between the fixed 
point of the trimer and those of the dimer which will turn out to persist also at the 
quantum level, allowing us to use the results reported in section |21 as a guide to gain 
deeper insight in the spectrum of the trimer. 

3.1. The mean- field picture of the trimer and its analogy to the symmetric dimer 

The effective Hamiltonian ensuing from Q,H t = U t 52 . I z j 1 4 + \w t ( | z\ \ 2 — \ z 2 \ 2 + \ z 3 | 2 ) — 
^T t [z 2 (zi + z 3 ) + z%{z\ + z£)], yields the following dynamic equations 

\hz 3 = (2U\z,\ 2 + ^) Zj - |z 2 , ihz 2 = (2U\z 2 \ 2 - ^) z 2 - I ( Zl + z 3 ) (7) 

where the first equation applies to both of the lateral sites, j = 1,3. This feature is a 

consequence of the mirror symmetry that Tit inherited from H t , and allows room for an 

integrable sub regime in the generally non-integrable dynamics of the trimer [T| I%1 130j . 

Notice indeed that the initial condition z\ = z 3 is propagated in time by (J7|). As we 

will discuss in some detail in the following, this sub regime of the trimer can be mapped 

exactly onto the asymmetric dimer considered in the previous section, which is why it is 

also referred to as dimeric sub regime of the trimer. However, since the analogy with the 

dimer applies more in general to the energies of most of the fixed points of the trimer, 

either exactly or approximately, we will refrain from assigning the adjective "dimeric" 

exclusively to the integrable sub regime of the trimer and the relevant fixed points. To 

avoid confusion, the latter will be referred to as symmetric fixed points (SFPs) alluding 

to their feature z\ — Z3. Likewise we will term asymmetric fixed point (AFP) any fixed 

point such that Z \ 7^ z 3 , with the only exception of the so-called central depleted well [T] 

(zi = — Z3, z 2 = 0) which, having no counterpart in the asymmetric dimer, will not be 

considered in the following discussion. 

The close relation between the integrable sub regime of the trimer and the 

asymmetric dimer can be emphasized by recasting equations (jZj), subject to the 

(dynamically propagated) constraint Z\ = z 3 , in terms of the alternative variables * 

d = (\zi\ 2 /2 — \z 2 \ 2 )/p t and (f> = (<px — 2 )/2.The resulting equations are exactly in the 

same form as those ensuing describing the dynamics of an asymmetric dimer with the 

following parameters 

3 1 
U d = -U t , w d = w t - -p t U t , T d = V2T t , p d = p t (8) 

Hence, as we mentioned in section the integrable sub regime of a trimer with no 
energy offset (wt = 0) could provide a realization of a rather asymmetric dimer, 

* Likewise, if the identification (zi 7 z r ) — (zi/\/2, Z2) is made, the resulting equations have the same 
form as those describing a dimer with site-dependent on-site interaction 1301 • 
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Figure 2. Trimer-dimer analogy for r t = 0.2. Left panel: exact mapping (JSJ) 
between the fixed points of the asymmetric dimer and the SFPs of the trimer. The 
gray symbols denote the local and stability character of the fixed point as evaluated 
in the entire phase space. More precisely gray stars, upward and downward triangles 
denote saddles, (either local or absolute) maxima and minima, respectively. White 
filled circles, dashed and solid thick black lines denote the same local and stability 
characters as evaluated within the sub manifold of the phase space pertaining to the 
integrable sub regime {z\ = 23). Right panel: approximate mapping © between the 
fixed point of the dimer and the AFP of the trimer. The gray squares denote the stable 
saddles typical of the non- integrable regime of the trimer whereas the remaining 
symbols are the same as in the left panel. The dotted black lines appearing in both 
panels signal the situation corresponding to a symmetric dimer. 

the corresponding relative offset being proportional to the (possibly very large) total 
population, u> d = —\ptU t . 

Here we will not give further analytic details about the study of the trimer dynamics 
(solution of the "asymmetric" fixed point equations, stability analysis), which can be 
found elsewhere |2] , our present concern being the analogy between the trimer and 
the dimer. In this respect we remark that the above described exact mapping applies also 
to the stability character of the SFP only within the "symmetric" phase subspace Z\ = z 3 
pertaining to the integrable sub regime. Indeed, when considering the entire phase space, 
the SFP corresponding to one of the (stable) maxima of the asymmetric dimer may turn 
out to be unstable. This is clearly shown in the left panel of figure EJ where we plotted 
the energies the SFPs, as obtained with the straightforward modifications to equation 
© ensuing from the mapping ©: u t (d;r t ) = ~ + |z/ d (<i;r d ), E t (d;r t ) = §£ d (d;r d ), 

As a further, somewhat unexpected, aspect of the dimer-trimer analogy, the 
parametric representation (JHJ) allows to describe also the energies of most of the 
asymmetric fixed points of the trimer, though approximately and for not exceedingly 



Persistence of mean-field features in the energy spectrum of small BEC arrays 9 



large values of r t . An example of this situation is portrayed in the right panel of figure 
121 where we plotted the energy of the AFPs (gray markers) along with those of the fixed 
points of an asymmetric dimer (black lines, circles) with exactly the same parameters 
as the trimer: 

U d = U t , w d = w t , T d = T t , Pd = Pt (9) 

As we mentioned above, the larger the rescaled hopping amplitude r t the less satisfactory 
the overlap between the fixed points of dimer and the AFPs of the trimer. However, the 
comparison among situations relevant to different hopping amplitudes seems to suggest 
that a more complex parameter mapping than © could produce better overlaps. Due to 
the limited space, here we do not develop this issue, but we rather limit our qualitative 
discussion to sufficiently small values of the hopping amplitude (up to ~ 0.5) for the 
simple mapping in Q to give satisfactory results. In this respect we mention that a closer 
analysis of the relevant fixed point equations j2] shows that the dimer analogy is more 
effective for the AFPs characterized by a very low population of one of the lateral sites. 
A further feature worth noticing is the correspondence between the dimeric minimum 
(thick solid line) and the stable saddle points Jj typical of the non-integrable regime of 
the trimer (gray squares). We conclude this somewhat qualitative analysis of the dimer- 
trimer correspondence by emphasizing that one of the dimeric maxima (black dashed 
lines) has no counterpart among the AFPs. As we will briefly discuss in the following, a 
similar feature appears in the in the purely quantum treatment of the problem, in terms 
of parity of the energy levels. 

3.2. Quantum persistence of the dimer-trimer analogy 

As we discussed at the beginning of the present section, the Hilbert space for a trimer is 
so large to require a numerical diagonalization of the (number conserving) Hamiltonian 
(JHJ) even for relatively small boson populations p t . Likewise, the pattern produced by 
the energy levels of the trimer as the relative energy offset w t is varied is much more 
complex than its dimeric counterpart (see e.g. figure EJ)- Some help in the classification 
of these energy levels comes from the observation that parity is a good quantum number, 
which casts some light on the rather complex avoided crossing pattern occurring as w t 
varies. The crossing of two levels is either forbidden or allowed depending on whether 
they have the same or different parity. Indeed levels with different parity belong to 
independent sectors of the Hilbert space, and there is no hindrance to their crossing. 

The discussion presented in section 13.11 provides the very interesting suggestion of 
interpreting the energy level pattern of the trimer in terms of the spectrum of a dimer 
with the same population and suitably chosen parameters. Actually, the dimer-trimer 
analogy persists also at the quantum level, as it is clearly shown in figure HH where we 
plotted the 36 levels of a trimer containing p t = 7 bosons along with the spectrum of two 

(J Recall that, according to the (linear) stability theory, the maxima and the minima are always stable, 
whereas the saddle point can be either stable or unstable. In the case of the trimer, the analytically 
predicted stability of the saddle points has been widely confirmed by numerical simulations [5]. 
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Figure 3. Quantum counterpart of figure [21 for r t = O.f and p t = 7. The thin 
solid lines and the circles denote even and odd energy levels of the trimer, respectively. 
The thick gray lines are the levels of the dimeric Hamiltonian . Those in the left 
panel have been obtained with the parameter choice © providing the exact mapping 
between the dynamics of Tld and the integrable sub regime of TLt- The dimeric levels 
in the right panel have been obtained according to the approximate mapping In- 
different dimers with the same population. The dimeric spectra appearing in the left 
and right panels of the figure were obtained from the diagonalization of Hamiltonian (JTJ) 
with the parameter choices in (jSJ) and (jHJ), respectively. Notice that in both cases each 
dimeric level follows quite closely a level of the trimer. In this respect we emphasize that 
the recognition of dimeric structures in the spectrum of © is quite difficult even after 
the introduction of an alternative canonical scheme incorporating the parity symmetry 
of the trimer [3U|, and involves non-trivial manipulations hardly conceivable except as 
an afterthought following the MF analysis. 

Notice that any definite-parity quantum state yields —n 3 ) = (afa^ — a^ai) = 0, 
which is a possible quantum counterpart of the MF integrability condition z\ = z^\\. 
This is basically the reason why each dimeric level in the right panel of figure El overlaps 
with an odd-even doublet of the trimer: the latter is the minimal set of eigenstates 
producing a state lacking a definite parity, which is the quantum counterpart of a MF 
asymmetric configuration. Recall indeed that the dimer under concern is the quantum 
counterpart of that qualitatively describing the asymmetric fixed points of the trimer. 
Conversely, though intrinsically asymmetric, the dimer considered in the left panel of 
figure El pertains to the symmetric configurations of the trimer and hence there is no 
surprise that in some instances a dimeric level corresponds to a solitary level of the 
trimer rather than to an even-odd doublet. If proved to occur in general, the fact that 

ft (a+a c ) = 0, where a c = {a\ + a 3 )/^/2 makes probably a better candidate. Due to limited space, we 
refer to jSOJ for a discussion about this issue. 
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such solitary levels are always even, as in figure remains to be understood. 

Notice that the (even) absolute maximum of the trimer for u t < seems to have 
a dimeric counterpart both in the left and in the right panel of figure El Actually, as 
it becomes evident considering a larger value of r t , the overlap in the right panel is 
accidental, and hence it is not an exception to our discussion above. This means that, 
as we mentioned in section |3~T1 one of the dimeric levels in the right panel (the absolute 
maximum for u t < 0) has no correspondent among the levels of the trimer. 

Owing to the quantum persistence of the dimer-trimer correspondence, most of 
the results about the asymmetric trimer presented in section El can be extended in a 
straightforward manner to the case of the trimer. 

As we mentioned at the beginning of the present section, figure |3] provides a quite 
useful key for interpreting the otherwise quite obscure trimeric version of figure 
Actually, owing to the quantum persistence of the dimer-trimer correspondence, most of 
the results about the asymmetric trimer presented in section |21 are so straightforwardly 
extended to the case of the trimer, that there is no need of including such figure in the 
present work. 

4. Conclusions and perspectives 

We focused on the Bose-Hubbard Hamiltonians describing two small arrays of interacting 
BECs, the asymmetric dimer and the trimer. In section |2] we related some interesting 
features of the energy spectrum of the former Hamiltonian, whose asymmetry ensues 
from an imbalance between the on-site energy terms, to the energies of the fixed point 
of the corresponding mean-field effective Hamiltonian. These results are inspired by and 
generalize those concerning the symmetric dimer reported in reference jTJ] . Section |3] is 
devoted to a similar investigation on the spectrum of the trimer. In order to find our 
bearings in the rather complex level pattern, in section 13.11 we exploited an interesting 
mapping among the mean-field dynamics of the trimer and that of the asymmetric dimer 
discussed in the previous section. In particular we showed that such mapping, which 
is exact only for the integrable sub regime of the trimer, applies approximately also to 
the energies of most of its fixed points. Based on these mean-field results, in section 
13.21 we observed that a similar approximate mapping exists also at the quantum level. 
We emphasize that the (quasi) dimeric character of the spectrum of the trimer is by 
no means evident from the structure of the relevant Bose-Hubbard Hamiltonian, and 
it strongly relies on the observations reported in section 13.11 The results gathered in 
sections 13.11 and 13 .2| along with the underlying systematic study of references [l] and 
[2], allow to extend straightforwardly the observations reported in section |21 to the case 
of the trimer. 

The results presented in this work form the basis for a systematic study of the 
purely quantum trimer. Among the issues worth developing we mention a possible 
characterization of the quantum counterpart of chaos. Recall indeed that, owing to 
the non-integrability of the mean-field dynamics, the trajectories based in the vicinity 
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of unstable fixed points are likely to display strong instabilities. Hence the mapping 
between fixed points and energy spectrum may provide a guide in recognizing possible 
quantum manifestations of chaos. 

Furthermore, if the system under concern is about to be realized experimentally, 
as promised by the recent advances in trapping techniques, the systematic study 
here proposed could provide a benchmark for verifying the effectiveness of the model 
Hamiltonians used to describe its physics. In this respect we recall that small bosonic 
lattices have been proposed for the implementation of quantum computation jHH IH2] . 
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